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Abstract 

Inertial sensors are commonly used to measure human head motion. Some sensors have been tested with dummy 
or cadaver experiments with mixed results, and methods to evaluate sensors in vivo are lacking. Here we present 
an in vivo method using high speed video to test teeth-mounted (mouthguard), soft tissue-mounted (skin patch), and 
headgear-mounted (skull cap) sensors during 6-13g sagittal soccer head impacts. Sensor coupling to the skull was 
quantified by displacement from an ear-canal reference. Mouthguard displacements were within video measurement 
error (<lmm), while the skin patch and skull cap displaced up to 4mm and 13mm from the ear-canal reference, 
respectively. We used the mouthguard, which had the least displacement from skull , as the reference to assess 6- 
degree-of-freedom skin patch and skull cap measurements. Linear and rotational acceleration magnitudes were over¬ 
predicted by both the skin patch (with 120% NRMS error for Umag, 290% for a^ag) and the skull cap (320% NRMS 
error for Omag, 500% for amag)- Such over-predictions were largely due to out-of-plane motion. To model sensor 
error, we found that in-plane skin patch acceleration peaks in the anterior-posterior direction could be modeled by an 
underdamped viscoelastic system. In summary, the mouthguard showed tighter skull coupling than the other sensor 
mounting approaches. Furthermore, the in vivo methods presented are valuable for investigating skull acceleration 
sensor technologies. 

Keywords: instrumented mouthguard, instrumented skin patch, instrumented skull cap, high speed video, soft tissue 
modeling 


1. Introduction 


and measurement accuracy. 


Traumatic brain injury biomechanics can be studied 
in human subjects using wearable head impact sensors 
that measure skull accelerations. The availability of 
low-power, low-cost MEMS accelerometers and gyro¬ 
scopes has spawned a flurry of head impact sensing 
approaches both in research and for consumer use. A 
helmet-mounted sensor system, the head impact teleme¬ 
try system (HITS), is an example of a widely used sens- 


ing option ( 

Duma et al. 

2005 1 Rowson et al. 2011 

[Rowson anc 

Duma[ 2011 

1 . However, factors such as 


helmet fit and padding type may affect sensor coupling 
to a human head, and in turn cause measurement errors 


([Higgins et al. 2007 1 Beckwith et al.||2012 Jadischke 


et al.j 201 3| l. More recently, industry and academic labs 
have developed alternative sensors with other form fac¬ 
tors and mounting locations such as the teeth, ear-canal, 
skin, and various types of headgear. For these devices, 
factors including fit, adhesion, soft-tissue elasticity, and 
hair/scalp properties may affect sensor skull coupling 


Instrumented bite blocks have been used in vivo as 
reference sensors ( [Funk et ak] [2009 1 jKnoxj |2004] l. A 
similar approach is to instrument a mouthguard, which 
is practical for held use in contact sports. Instru¬ 
mented mouthguards have been evaluated in vitro with 


a clamped-jaw dummy as a reference (Camarillo et al. 


|2013] l. One error source was introduced when the sen¬ 
sors were placed in a protruding tab on the mouthguard, 
which exhibited a mechanical resonance. The resonance 
led to errors in peak acceleration measurements, but 
RMS acceleration errors were still within 10%. |Bartschl 


et al. (20141 also demonstrated accuracy using a dummy 


head that does not have a lower jaw to clamp the mouth¬ 
guard. 


Skin patch and skull cap sensors are also used in re¬ 
search and are becoming commercially available. Skull 
coupling of either approach has not been evaluated in 
literature. Previous studies with skin-mounted sen¬ 
sors and optical markers at other locations on the body 
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(e.g. knee joint) report artifacts due to skin dynamics 
( |Reinschmidt et ari|1997[[Lucchetti et al.[|1998HShult^ 


|et al. 201 1| |. Differentiation of position and orienta 
tion measurements can amplify acceleration errors from 
such soft tissue artifacts. Therefore, some researchers 
have modeled soft-tissue dynamics to correct for mea¬ 
surement errors ([Trujillo and Busb^ 1990| Kim et al. 


|1993) . It is likely that soft tissue or textile dynamics af¬ 
fects sensor performance, but the effects have not been 
quantified for head impacts. 

Sensors are commonly evaluated through in vitro (an¬ 
thropomorphic test device) and/or ex vivo (postmortem 
human subject) methods. Both methods have the ad¬ 
vantage of using high fidelity reference sensors rigidly 
attached to the skull, methods test sensor accuracy with¬ 
out confounding factors such as soft tissue motion, and 
are valuable for verifying sensor hardware selection and 
programming. Ex vivo methods introduce additional 
biofidelic factors including skull and tissue dynamics, 
but the effects of postmortem changes in tissue proper¬ 
ties are unknown, and the lack of muscle forces may 
also affect head dynamics. Thus in vivo sensor evalua¬ 
tion may help to account for these factors. But there is a 
lack of in vivo methods, since we can neither screw ref¬ 
erence sensors directly to the skull, nor use dangerous 
impact conditions in human subject studies. 

Thus our objective is to develop a non-invasive in vivo 
method to evaluate head impact sensing approaches. 
Using this method, we will test three types of ap¬ 
proaches: sensors fit to hard tissue (teeth), adhered to 
soft tissue (skin), or mounted on headgear that fits the 
head. Corresponding to these sensor types, we assessed 
skull coupling of instrumented mouthguard, skin patch, 
and skull cap sensors in a human subject during mild 
soccer head impacts using high speed video tracking 
(Fig. [T] Methods Sec 2.3, Results Sec 3.1). Then we 
used the sensor with the least amount of skull displace¬ 
ment, which is the mouthguard, as the reference sen¬ 
sor to compare full 6-degree-of-freedom (6DOF) sen¬ 
sor measurements (Methods Sec 2.4, Results Sec 3.2). 
To account for sensor error from soft-tissue motion, 
we modeled sensor-skull dynamics using a simple vis¬ 
coelastic model (Methods Sec 2.5, Results Sec 3.3). 


2. Methods 

2.7. Human Subject Experiment Setup 

A 26 year-old male human subject underwent soccer 
head impacts with clenched teeth, at initial ball speed 
of 7m/s (Fig. [^), which is the average header speed 
in youth soccer, and higher than that (5.7m/s) in adults 


( Shewchenko et aT| [2005 | l. A ball launcher (Sports 
Tutor, Burbank, CA) helped to simulate a kicked ball, 
and the ball was inflated to approximately 8-9psi. Hu¬ 
man subject protocols in this study have been approved 
by the Stanford Institutional Review Board (IRB No. 
26620), and informed consent was obtained from the 
subject. 

2.2. Instrumentation 

volunteer wore a custom-fit mouthguard ( jWuj 
2014| l, a skin patch adhered to skin on the mas¬ 
toid process (xPatch Gen2, X2Biosystems, Inc.), and an 
elastic skull cap (Reebok). The mouthguard had ap¬ 
proximately 4mm average thickness and 7mm height 
above the gum line. Electronics were placed inside 
the mouth to avoid tab resonance, while video mark¬ 
ers were fixed on a light-weight protruding tab. The 
volunteer’s head circumference measured 60.3cm, and 
wore a size large elastic skull cap. The same electron¬ 
ics from the mouthguard were placed on a soft card¬ 
board in the lateral insert of the skull cap. To enable 
comparison of measurements, devices were set to a low 
triggering threshold (4 - 6g) to ensure all impacts were 
recorded on each device. Sensor signals were synchro¬ 
nized through videoby aligning sensor-measured kine¬ 
matics with video-derived sensor kinematics for each 
individual sensor. To derive sensor kinematics from 
video, we resolved the position and orientation of each 
sensor using a previously-described method of tracking 
fiducial markers in stereo video |Hernandez et al.| ( |2015l l. 
During synchronization, we found that skin patch mo¬ 
tion can lag behind the mouthguard by as much as 15 
ms. So we recorded 30ms of pre-trigger data and 70ms 
post-trigger. 

In addition to these three sensors, the subject wore 
a deeply-inserted earplug as a skull reference point. 
Previous research confirmed skull coupling of deeply- 
inserted custom-fit silicone earplugs using reference 
sensors screwed onto the skull of a postmortem human 
subject ( jSalzar et al.||2014l|Christopher et al.||2013| l. An 
expandable foam earplug was inserted approximately 
20mm into the ear canal, similar to the depth at which an 
ear sensor in the postmortem experiment was mounted. 
The low mass of a foam earplug (0.2g, compared to 5g 
for a typical custom-formed silicone earplug) also min¬ 
imizes inertial effects and improves coupling to the ear 
canal. This skull reference was used to assess skull cou¬ 
pling of the three sensors. 

2.3. Video analysis of skull coupling 

We took high speed stereo video at 1000 frames per 
second and 1920x1200 resolution (0.3mm/pixel at dis- 
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tance of head), using two Phantom Miro LC320 cam¬ 
eras (Vision Research, Wayne, NJ), to track motion of 
the sensors (Fig. By Fiducial markers with 1.5mm 
square grids were fixed onto each sensor and the earplug 
to track the change in distance between grid centroids. 
The mouthguard, skin patch, and skull cap sensors each 
had 26-30 trackable points, and the earplug had 1-3 
trackable points. The two cameras were positioned such 
that 1) the tracking grid on the deeply inserted earplug 
had at least 1 trackable point and 2) all three sensors 
were visible throughout the head impact. Due to these 
constraints, there was a triangulation angle of 7.4 de¬ 
grees between the cameras. Using the Camera Calibra¬ 
tion Toolbox for Matlab (Bouguet 2013t Zhang 1999 


Heikkila and Silven 1997| l, we performed stereo call 

bration to enable 3D position tracking in the head mo¬ 
tion space. 

To verify our video method, we tracked the 4 corners 
of a 20cmx20cm calibration grid during 6DOF motion 
through 4000 frames (4 seconds). The calibration grid 
moved through the same region of space as head motion, 
at the same distance away from the cameras as the sub¬ 
ject. (R3-4) We assessed both planar and depth measure¬ 
ments of distance, by comparing stereo video measure¬ 
ments with ground truth grid distances in each frame. In 
addition, we derived sagittal kinematics of each sensor 
from video measurements, to cross validate with sensor 
measurements. Position time histories of points fixed on 
the grid were combined to determine the least-squares 
rotation matrix describing the orientation of the body- 
fixed frame of each sensor in the camera-fixed frame 
( [Hernandez et al.|[2015| l. The time derivative of this ro¬ 
tation matrix was related to the rotational velocity of 
the body-fixed frame in the camera-fixed frame. We 
used this rotational velocity to take the time deriva¬ 
tive of linear position and velocity to determine veloc¬ 
ity and acceleration, respectively, in the moving body- 
fixed frame. Due to near-parallel arrangement of cam¬ 
eras, out-of-plane measurements are expected to have 
larger errors, and we compared 3DOF sagittal kinemat¬ 
ics instead of full 6DOF kinematics when validating our 
video method. 


2.4. Comparison of6-DOF sensor kinematics 

Video measurements of the ear-canal reference only 
allowed for 1-3 trackable points, so 6DOF kinematics 
could not be computed. To evaluate 6DOF measure¬ 
ment differences, we selected a reference device with 
the least relative motion to the ear-canal - the mouth¬ 
guard sensor. 

We transformed kinematic measurements of the skin 
patch and skull cap sensors to an estimated center-of- 


gravity (CG) location, for trials where mouthguard and 
skull were found to have the lowest relative displace¬ 
ment (see Results Sec. 3.1). CG location was estimated 
based on a 50th percentile male human head model. We 
first estimated the projection vector from mouthguard 
to CG by using the upper dentition as an anatomical 
landmark on the model. Then we measured the pro¬ 
jection vectors of the skin patch and skull cap sensors 
to the mouthguard location on our human subject. Us¬ 
ing stereo video of the subject standing still without 
head motion, we derived relative position and orienta¬ 
tion of the 3 devices and confirmed with physical mea¬ 
surements. The resultant projection vectors for the skin 
patch and skull cap sensors were calculated by summing 
their projections to mouthguard and the projection from 
mouthguard to CG. Using this method, we ensure that 
all three devices can be projected to the same point in 
space for comparison, even though the model-estimated 
CG location is expected to have some error. 

Linear accelerations were projected to the CG loca¬ 
tion for comparison using the following rigid body vec¬ 
tor relationship; 

dcG - Os + d y. Ts + o) X (o) X r^) (1) 

where acc is head linear acceleration at CG, a^ is head 
linear acceleration at each sensor location, a is head an¬ 
gular acceleration measured by the sensor, u is head 
angular velocity measured by the sensor, and r^ is the 
vector position of CG location from the sensor location. 


Transformed skin patch and skull cap sensor data 
were compared with mouthguard reference data in 
the anterior-posterior (AP), left-right (LR), inferior- 
superior (IS) directions for linear acceleration, and the 
coronal, sagittal, and horizontal planes for rotational ac¬ 
celeration. Quantities reported include vector magni¬ 
tudes and individual axis differences from all 6DOF lin¬ 
ear acceleration, angular velocity, and angular acceler¬ 
ation. We performed linear regression analysis of peak 
kinematic values. We also reported the average devi¬ 
ation in peak values from the mouthguard reference, 
and compared the directions of head acceleration. In 
addition to peak values, we assessed the agreement of 
time traces by computing root-mean-square (RMS) dif¬ 
ference and normalized root-mean-square (NRMS) dif¬ 


ference for 25 samples around peak measurements (Ca- 


jmarillo et al.j [201 3| l. To better understand the sources 
of sensor errors prior to projection to CG, we also com¬ 
puted the RMS and NRMS differences at the sensor lo¬ 
cation. That is, we compared skin patch/skull cap sig¬ 
nals at the skin patch/skull cap location with mouth- 
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guard reference signals projected to the skin patch/skull 
cap location. 

2.5. Modeling of sensor dynamics 

2.5.1. Model description 

Since the soccer headers were frontal hits and head 
motion was mainly anterior-posterior, we modeled the 
skin patch and skull cap translation in the anterior- 
posterior direction as a second order linear system with 
base (skull) excitation (Fig. |^, Equation]^, where 

nisensor ' Combined mass of the sensor and soft tis¬ 
sue (kg) 

Kt - effective linear spring constant of sensor-skull 
mounting (N/m) 

Ct - effective linear damping constant of sensor- 
skull mounting (N-s/m) 

dsensor “ absolute displacement of the sensor (m) 
dskull - absolute displacement of the skull (m) 

sensord sensor — ^t(dsensor dskull) Ctidsensor dskullli'^') 

We also modeled sagittal rotation of the head using a 
similar second order linear system (Fig. |2f5, Equation 

0 - 

^sensor^sensor ~ sensor' -Sskull)-Cr(0 sensor' -OskuldO) 

2.5.2. Model evaluation 

We used sensor data from soccer impacts to hnd 
the input-output relationship defining the skull-tissue- 
sensor system. Only trials where mouthguard-skull cou¬ 
pling was best (<0.5mm) were modeled. We ht model 
parameters using mouthguard measurements as skull in¬ 
put, and skin patch/skull cap measurements as sensor 
output. 

Sensor modeling was done by htting spring, damper, 
and mass parameters to the input skull kinematics and 
output sensor kinematics. This system takes prescribed 
skull motion (obtained by transforming the mouthguard 
kinematic signal to the sensor location) as an input and 
does a forward dynamics simulation using the Matlab 
ode45 integrator to solve for the resulting sensor kine¬ 
matics. The initial state for the system is set such that 
there is no relative motion between the skull and sensor 
and the spring is at its rest length. We used the Matlab 
function fmincon to find a set of spring stiffness, damp¬ 
ing coefficient, and mass parameters that minimize the 
root mean squared error between the output model sen¬ 
sor kinematics and the measured sensor kinematics. The 


search space for these parameters are bounded to con¬ 
strain the system to realistic values and the initial guess 
is set as the midway point between bounds. Considering 
that the sensors moved along with packaging and under¬ 
lying soft tissue, a loose mass bound of 1 to 50 grams 
was placed on the optimization to determine the mass of 
the system. Each impact is separately fitted in order to 
account for variabilities in impact conditions (location 
and force). The NRMS of the fit was used to assess the 
model. 

In order to check the validity of the linearity assump¬ 
tion in our model, we examined the input-output rela¬ 
tionship in the frequency domain by plotting the ex¬ 
perimental and theoretical frequency response functions 
(FRFs). Experimental FRFs were calculated by taking 
the ratios of the Fast Fourier Transform (EFT) ampli¬ 
tudes of the skull input and skin patch output. Fur¬ 
thermore, the analytical frequency response of the base- 
excitation model depicted in Fig. was compared 
against the experimental FRFs. 


3. Results 

3.1. Skull coupling from video 

We verified stereo video tracking to have <lmm er¬ 
ror in the sagittal plane (Fig. [^). When the calibration 
grid was displaced or rotated in a plane perpendicular to 
the camera axis (sagittal), stereo triangulation estimated 
distances between points on the calibration grid with 
<lmm error. When depth measurement was involved, 
such as during out-of-plane rotation (i.e. grid was tilted 
from camera plane), errors were slightly larger but still 
within 2mm. Overall, more than 90% of the errors were 
within 0.5mm. As an additional verification step, we 
compared video-derived sagittal kinematics of each sen¬ 
sor with those measured by the accelerometer and gyro¬ 
scope in each sensor. These kinematics matched with 
<30% NRMS difference, even for linear acceleration 
double-differentiated from position tracking (Fig. Up, 
Table[T}. 

For the mouthguard, all head impact trials (n=16) 
showed relative displacements from the earplug of 
<lmm (Fig. p=0.5mm, cr=0.2mm), within video 
measurement error. In addition, 10 of the 16 trials 
showed relative mouthguard displacements of within 
0.5mm. In contrast, the skin patch sensor displaced 
by 2-4mm (//=3mm, cr=0.7mm) at the moment of head 
impact; the skull cap sensor displaced by 2-13mm 
(//=5mm, cr=3mm). 
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3.2. 6-DOF sensor kinematics 

Using the mouthguard as a reference, head accelera¬ 
tion peak magnitudes for the 10 trials averaged 9.3+2g 
in linear acceleration and 75Q±3QQradI s^ for rotational 
acceleration. Skin patch estimation of head linear and 
angular acceleration peak values were over-predicted by 
15+7g and 2500+1200ra£//i^ on average, respectively; 
the skull cap values were over-predicted by 50+3 Ig and 
4300±2700radls^ (Tablej^. These over-predictions are 
the average differences between the sensor measure¬ 
ments and the reference (mouthguard) measurements 
Fig. 1^ also shows the over-predictions in peak mag¬ 
nitudes, with patch/cap predictions scattered above and 
away from the m=l reference line. In addition, patch 
and cap peak vector magnitudes and all individual com¬ 
ponent peaks had large variances, and did not cotTelate 
well with the mouthguard reference (Supplemental Fig. 
IST) . When y-intercept is forced to be 0 in the regression, 
the coefhcient of determination is sometimes negative 
(Supplemental Table |Sl| l, indicating that the residual er¬ 
ror in the fit is greater than the variance in the data. To 
evaluate dynamic relationships among sensors. Table 
reports the time lag/lead between the mouthguard peak 
value and patch/cap peak values. The skin patch linear 
acceleration magnitude peak was the most consistent, 
occurring 15±3ms after the mouthguard peak. 

The skin patch and skull cap also predict different 
directions of head motion compared with the mouth¬ 
guard. Fig. 1^ and Table|^show the differences in kine¬ 
matic vectors at the moment of peak magnitude. The 
mouthguard reference measured head motion to exhibit 
mostly planar motion with anterior-posterior (AP) lin¬ 
ear acceleration and sagittal rotation. In contrast, skin 
patch motion was mostly out-of-plane with left transla¬ 
tion and horizontal rotation, and the skull cap was not 
in a consistent direction. Breaking magnitudes down 
into per-component comparisons. Fig. shows sam¬ 
ple 6DOF kinematic waveforms of a representative im¬ 
pact with mostly anterior-posterior motion. We observe 
over-predictions of kinematics in all axes, including out- 
of-plane (non-sagittal) directions. In fact, the highest 
peaks for the skin patch all occur in out-of-plane axes: 
left-right (LR) linear acceleration and horizontal rota¬ 
tion. The skin patch signals also show damped oscilla¬ 
tory behavior. Table shows RMS and NRMS errors 
of the skin patch and skull cap in 6DOF. All magnitude 
errors were above 100%. At sensor location, the lin¬ 
ear acceleration RMS and NRMS errors are lower than 
those at the CG (comparing Table and Supplemen¬ 
tal Table S2i. Linear acceleration magnitude has 18% 
NRMS error for the skin patch, and 60% NRMS error 
for the skull cap, compared to 120% and 320% when 


projected to CG. Angular velocity and angular acceler¬ 
ation are independent of the location on the head, and 
have the same high errors at both sensor location and 
CG. . 


3.3. Sensor modeling 


Identified model parameters are detailed in Table 
including the NRMS errors of the fit. The statistics are 
generated from: all trials for skin patch AP translation; 
7 of 10 trials for skin patch sagittal rotation; and 8 of 
10 trials for cap AP translation. Since this is a very 
simplified model, we did not expect all trials to fit to 
the model. In order to report meaningful parameter esti¬ 
mates, we omitted trials where the model could not fit to 
the data. Some trials were omitted for patch sagittal ro¬ 
tation since the input mouthguard signal lagged behind 
the patch measurement. In such cases, the mouthguard 
signal could not be used as input for modeling. Instead, 
the input to the system likely came from other coupled 
axes. Some trials were omitted for skull cap AP trans¬ 
lation, since reasonable estimates of model parameters 
could not be found due to sharp spikes in the signal, 
which likely resulted from direct soccer ball impact. For 
all trials of skull cap sagittal rotation, the mouthguard 
signal consistently lagged behind the skull cap signal, 
and the model failed to fit. Fig. shows sample traces 
comparing the model fit and measured signals. Over¬ 
all, model fits were the most consistent for skin patch in 
AP translation. The average estimated mass of the skin 
patch system (8.5g) is greater than mass of the sensor it¬ 
self (5g), showing that some underlying tissue mass was 
also activated during the impacts. 

For skin patch AP translation, where the model could 
be fit to all 10 trials with relatively low variance in 
model parameters, we verified the linear system as¬ 
sumption and examined the frequency response around 
the resonant frequency of the model. Fig. and B 
show the FFT amplitudes for the mouthguard (skull) in¬ 
put and skin patch (sensor) output, respectively. The 
FRFs, which are the ratio of output FFT amplitude to the 
input FFT amplitude, are shown for all trials in Fig. |^. 
We compared these with a theoretical FRF (red line), 
which is the frequency response of the transfer func¬ 
tion for the model in Fig. with the average model 
parameters (Table |^. As shown, 9 of the 10 trials had 
FRFs with amplitudes and peak frequencies within a +1 
standard deviation confidence interval, which validates 
the linearity assumption, since a linear system exhibits 
a consistent FRF ( |Ewins[^000| l. 
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4. Discussion 

In this study, we developed an in vivo method to 
quantify skull coupling of teeth-mounted (mouthguard), 
soft tissue-mounted (skin patch), and headgear-mounted 
(skull cap) sensors. All mouthguard-skull displace¬ 
ments were within video measurement error. This was 
an expected outcome for mild impacts with clenched 
teeth, since the mouthguard is custom-formed to the 
teeth. It does conhrm that in vivo clenching force can 
hold the mouthguard in place for the soccer impact con¬ 
ditions tested. On the field, other practical factors such 
as mandible motion and variations in mouthguard fabri¬ 
cation may affect performance. Unanticipated impacts 
with open jaw may also pose a problem for the mouth¬ 
guard. Thus held mouthguard coupling remains to be 
tested. However, in anticipated impacts, the teeth are 
likely clenched, resulting in similar conditions as our 
tested scenario. 

Non-rigid skull coupling led to skull measurement 
errors in skin patch and skull cap sensors. Both sen¬ 
sor magnitudes showed over-predictions of peak head 
CG kinematics. Such over-predictions resulted from 
measurements of signihcant out-of-plane motion, while 
head motion was mostly sagittal. In fact, the skin patch 
and skull cap often measured acceleration peaks in a dif¬ 
ferent vector direction from the mouthguard reference 
with >50 degrees deviation (Table |^. With the over¬ 
predictions and differences in direction, raw data from 
these sensors likely cannot be directly used to predict or 
study injury risks. Sensor errors need to be corrected via 
models, and/or reduced by improving skull coupling. 

In the primary plane of motion (sagittal), a sim¬ 
ple viscoelastic model approximated single-degree-of- 
freedom skin-sensor dynamics, in agreement with pre¬ 


vious studies of lower-limb soft tissue dynamics (Tru 


jillo and Busby 1990t Kim et al. 1993| l. As a feature 
of a linear dynamic system, the time lag between skin 
patch translation and skull translation was consistently 
around 13ms (Table |^. This time lag was due to vis¬ 
coelasticity of the skull-skin-sensor system. Skin patch 
motion in this axis had linear model hts with low vari¬ 
ance in parameters (Table ffl. From the frequency re¬ 
sponse of the system (Fig. we further conhrm the 
underlying linear dynamics of the system. For the soc¬ 
cer impact, the skull input excited a frequency range that 
includes the resonant frequency of the skull-skin-sensor 
system (20-30Hz), where the gain of the system is max¬ 
imum. However, if the head is driven at a different input 
frequency (i.e. a different input duration) and/or at con¬ 
siderably higher amplitudes, the peak gain may vary, or 
a different mode of the system may be excited. There¬ 


fore, it would be an oversimplihcation to use a static 
gain term to estimate head acceleration. 

Behavior of the skull cap sensor, especially in rota¬ 
tion, was less predictable than the skin patch sensor, 
varying from impact to impact. The time lag between 
cap and skull motion had large variance (Table|^, indi¬ 
cating inconsistent behavior. This is likely due to direct 
impact of the skull cap by the ball. Ball force and im¬ 
pact location were important factors for the skull cap 
system not accounted for by the simple linear model. 
During our experiments, the skull cap sometimes com¬ 
pletely dislocated from the head (results not included 
since skull cap was not in camera view). However, it 
is possible that for other impact conditions, such as hel- 
meted impacts, the skull cap may have less relative mo¬ 
tion from the skull. 

Although the modeling method shows promise in 
simulating soft tissue behavior, this may not be suffi¬ 
cient to correct for sensor error. Skin patch linear accel¬ 
eration in the anterior-posterior direction had the lowest 
variance in the optimized model parameters. To esti¬ 
mate skull input from skin patch measurements, an in¬ 
verse dynamics simulation can be performed using av¬ 
erage model parameters. But this would only help to 
mitigate skin patch errors in IDOF, and the variability 
in tissue response across different impact conditions and 
subjects can pose a challenge in developing a universal 
model. Also, referring to Supplemental Table S2 we 


show that the linear acceleration errors at head CG are 
much higher than those at the sensor locations. This 
indicates amplihcation of rotational velocity and accel¬ 
eration errors when measurements are projected to CG. 
Thus better models or design changes to mitigate rota¬ 
tional errors may help significantly improve sensor per¬ 
formance. 

The in vivo methods in this study have some limi¬ 
tations. First, only mild impacts were assessed in a 
single human subject. We tested low-speed 7m/s im¬ 
pacts for protection of the human subject, while held 
ball speeds could reach up to 17m/s ( |Shewchenko et al^ 
2005]l. Also, with one subject and one impact condi¬ 


tion, the outcomes are likely subject to variability in 
soft tissue properties, skull cap ht, mouthguard ht, and 
impact location/severity. The head acceleration levels 
(6-13g) are also low compared with injury-level accel¬ 
erations on the held ( [Hernandez et aH 2014j l, for pro¬ 
tection of human subjects. Second, high speed stereo 
video tracking was limited by the need to deeply insert 
the ear-canal reference for tight skull coupling. This 
led to near-parallel arrangement of cameras and low 
number of trackable ear-canal points. As a result, we 
could not derive 6DOF reference measurements from 


6 




















the skull reference. Third, we made the assumption 
that the ear canal skull reference is rigidly attached to 
the skull. This limitation is difficult to eliminate in an 
in vivo study. However, the fact that the mouthguard 
and the ear reference showed low relative motion gives 
confidence by cross validating these two mounting lo¬ 
cations approximately 10cm apart on the head. Fourth, 
mouthguard bite force was not controlled or measured 
in the experiment. At the low acceleration levels in this 
experiment, we do not expect bite force to signihcantly 
change results. However, at higher accelerations, bite 
force may need to be quantihed. Lastly, in practice, 
higher degree-of-freedom dynamic models with more 
elements/parameters may be necessary to predict both 
in-plane and out-of-plane sensor errors. 

In summary, we have developed a method to quantify 
skull coupling of wearable head impact sensors in vivo, 
and evaluated some common sensing approaches. The 
instrumented mouthguard was shown to have close skull 
coupling when clenched during mild soccer head im¬ 
pacts. The skin patch and skull cap devices had higher 
displacements from the skull. Raw data from sensors 
without close skull coupling should be interpreted cau¬ 
tiously both in trauma research and clinical assessment. 
To mitigate insufficient coupling, design modihcations 
and modeling may help to reconstruct skull motion. 
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Table 1: RMS errors of video-tracked sagittal kinematics with respect to sensor measurements 


AP linear acceleration IS linear acceleration Sagittal angular velocity 


Device 

RMS (g) 

NRMS (%) 

RMS (g) 

NRMS {%) 

RMS (rad/s) 

NRMS (%) 

Mouthguard 

1.0 (0.6) 

15.5 (6) 

1.8 (1) 

18.1 (10) 

0.6 (0.3) 

12.2 (7) 

Skin Patch 

1.4(0.3) 

13.7 (2) 

1.3(1) 

29.3 (30) 

2.1 ( 0.7) 

28.6 (9) 

Skull Cap 

4.4 (3) 

16.0 (9) 

2.1(1) 

12.5 (5) 

2.0 (2) 

13.2 (11) 


Note: Average RMS and NRMS errors are reported with standard deviation in parentheses. 
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Table 2: Comparing skin patch and skull cap kinematic peak values with mouthguard reference 


Linear Acceleration Magnitude 


Avg Peak 

Std Peak Diff 

Avg Peak 

std Peak 

Avg Direction 

Std Direction 

Device 

Diff (g) 

(g) 

Delay (ms) 

Delay (ms) 

Diff (deg) 

Diff (deg) 

Skin Patch 

+15 

±7 

+ 13 

± 2 

123 

±26 

Skull Cap 

+ 50 

±31 

+ 4 

±4 

54 

±38 

Angular Velocity Magnitude 


Avg Peak 

Std Peak Diff 

Avg Peak 

Std Peak 

Avg Direction 

Std Direction 

Device 

Diff (rad/s) 

(rad/s) 

Delay (ms) 

Delay (ms) 

Diff (deg) 

Diff (deg) 

Skin Patch 

+ 9.9 

±4 

-7 

± 15 

51 

± 11 

Skull Cap 

+ 10.3 

±8 

-9 

± 12 

58 

±22 

Angular Acceleration Magnitude 


Avg Peak 

Std Peak Diff 

Avg Peak 

Std Peak 

Avg Direction 

Std Direction 

Device 

Diff (rad/s^) 

(rad/s^) 

Delay (ms) 

Delay (ms) 

Diff (deg) 

Diff (deg) 

Skin Patch 

+ 2500 

± 1200 

+ 15 

±9 

116 

±22 

Skull Cap 

+ 4300 

± 2700 

+ 4 

±3 

98 

±27 


Note: Positive (+) sign indicates over-prediction for average peak difference, and indicates delay for 
average peak delay. 


10 











Table 3: Sensor RMS differences with respect to mouthguard time traces at head CG 




Linear Acceleration 

Angular Velocity 

Angular Acceleration 

Device 

Parameter 

RMS, g 

NRMS, % 

RMS, rad/s 

NRMS, % 

RMS, rad/s ^ 

NRMS, % 

Skin Patch 

AP Translation/ 
Coronal Rotation 

12.5 (4) 

120 (30) 

1.2 (0.6) 

78 (29) 

182 (114) 

44 (27) 


LR Translation/ 
Sagittal Rotation 

8.8 (3.3) 

300 (160) 

2.2 (0.6) 

48 (16) 

786 (170) 

120(62) 


IS Translation/ 
Horizontal Rotation 

5.1 (1.6) 

70 (30) 

7.5 (2.2) 

700 (380) 

1772 (516) 

860 (500) 


Magnitude 

10.8(4.5) 

120 (40) 

5.7 (1.9) 

150(76) 

1539 (552) 

290 (230) 

Skull Cap 

AP Translation/ 
Coronal Rotation 

12.4(6.8) 

140 (110) 

4.8 (2.5) 

370 (280) 

1858 (1290) 

460 (340) 


LR Translation/ 
Sagittal Rotation 

15.7(8.3) 

510 (300) 

4(2.5) 

83 (52) 

1349 (823) 

160(60) 


IS Translation/ 
Horizontal Rotation 

19 (13.2) 

250 (180) 

2.6 (1.5) 

330 (450) 

945 (565) 

590(660) 


Magnitude 

25.6 (16.1) 

320 (260) 

4.6 (3.5) 

110 (80) 

2233 (1379) 

500(140) 


Note: average RMS and NRMS errors are reported \A/ith standard deviation in parentheses. 
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Table 4: Model parameters 


AP Translation 

Device 

Kt (N/m) 

C t (N-s/m) 

^ sensor (Q) 

f(Hz) 


NRMS (%) 

Skin Patch 

189.4 (46) 

0.67 (0.16) 

8.5 (1.6) 

23.8(2.3) 

0.27 (0.09) 

20.4(11) 

Skull Cap 

396 (125) 

1.43 (0.40) 

12.4 (9.6) 

35.2 (16) 

0.40 (0.17) 

31.1 (11) 

Sagittal Rotation 

Device 

K, (N-m/rad) 

C ^ (N-m/rad-s) 

1 sensor ) 

f(Hz) 


NRMS (%) 

Skin Patch 

0.30 (0.08) 

4x 10 '‘(lx lO '^) 

4.5x 10® (lx 10®) 

41.6 (5) 

0.20 (0.08) 

20.0 (8) 

Skull Cap 

- 

- 

- 

- 

- 

- 
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Figure 1: In vivo evaluation and comparison of instrumented mouth¬ 
guard, skin patch, and skull cap. (A) A human subject underwent mild 
soccer head impacts, wearing all three sensors. Fiducial markers were 
mounted on the head, with one set on a deeply-inserted earplug (skull 
reference), and a set on each sensor. (B) Markers were tracked us¬ 
ing high-speed stereo video to determine the relative motion between 
each sensor and the skull reference. 
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Figure 2: Dynamic model description. (A) For anterior-posterior 
translation, we modeled each sensor-skull system as a second order 
linear system with a spring and a damper in parallel. These elements 
represent dynamics of the underlying tissue as well as the packaging 
and attachment of the sensors. (B) For the case of sagittal rotation, the 
elements are torsional springs and dampers. 
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Figure 3: Video validation. (A) We calculated errors from video tracking using a 20cm x 20cm calibration grid moving through the head motion 
region of space. When the calibration grid displaced or rotated with a sagittal orientation (i.e. planar measurements), errors were always sub¬ 
millimeter. When depth measurement was involved, with the grid rotating in non-sagittal directions, errors were larger but still within 2mm. (B) 
We fit a t location-scale distribution to the error, and there is less than 2.5% total probability of errors greater than 1mm. (C) We also verified that 
video-derived sagittal kinematics agree well with those measured by the sensors, which further confirms our video measurements. 
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Figure 4: Sensor coupling evaluation. Using high speed video, we compared the relative displacements the three sensors from the skull. Among 
16 trials, the mouthguard always had sub-millimeter displacements from the skull within video error, while the other two sensors had higher 
displacements. 
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Figure 5: Compaiing magnitude and direction of peak kinematics. Skin patch and skull cap measurements were compared with the mouthguai'd, 
which was used as the skull reference. (A) shows scatter plots of peak magnitudes of linear acceleration, angular velocity, and angular acceleration. 
The m=l reference line represents the ideal correlation when skin patch and skull cap match the reference. In (B), we plot vectors showing the 
direction of head accelerations and velocities at the moment of peak magnitude. The mouthguard, or skull reference, measured head motion 
to exhibit mainly posterior or superior linear acceleration, with sagittal rotation. The other two sensors, however, predict different directions of 
acceleration/velocity. 
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Figure 6: Comparing 6DOF kinematics in a sample impact. Measurements from skin patch and skull cap are compared with the mouthguard 
reference. Both sensors over-predict accelerations in the sagittal plane (highlighted axes) as well as out-of-plane axes. 
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Figure 7: Model predictions for skin patch and skull cap. For skin 
patch AP translation, sagittal rotation, and skull cap AP translation, 
we could fit underdamped second order linear systems to model sensor 
output. For skull cap y rotation, the mouthguard (skull) input lagged 
behind the sensor output, and thus this axis could not be modeled 
using this simple model. 


19 
















10 50 100 150 200 



10 50 100 150 200 

Frequency (Hz) 



Frequency Response Function 


Empirical trials 



Figure 8: Frequency response function of skin patch in AP linear 
acceleration. The FFT of the mouthguard input (A) and skin patch 
output (B) both show peak amplitudes occumng in a low frequency 
range. The frequency response function of the system (C) shows that 
for 9 of 10 of the trials modeled, the frequency response functions are 
similar, within 1 standard deviation (shaded region) of the theoretical 
model (red line), which further demonstrates lineality of the system. 
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Table SI: Linear Regression of Peak Kinematics 




Linear Acceleration 

Angular Velocity 

Angular Acceleration 

Device 

Parameter 

m 

1 

r 

m 

1 

r 

m 

2 

r 

Skin Patch 

X 

-0.7 

0.01 

0.8 

-0.09 

0.2 

0.00 


y 

-3.0 

-9.40 

1.5 

0.50 

-0.9 

-0.29 


Z 

0.5 

0.06 

5.3 

-10.11 

-0.9 

-4.60 


magnitude 

2.6 

0.15 

3.1 

-1.02 

3.5 

-3.41 

Skull Cap 

X 

3.6 

-0.78 

5.7 

-1.07 

-9.0 

-0.08 


y 

-2.3 

-1.60 

2.2 

0.14 

-0.8 

0.03 


Z 

6.1 

0.17 

2.4 

0.30 

-0.8 

-2.06 


magnitude 

6.0 

-0.28 

3.4 

0.19 

6.7 

0.56 


* Linear regressions are forced to have zero y-intercept. A negative coefficient of determination (r^ ) indicates that the fit is worse than a 
horizontai iine passing through the mean, and a constant term is needed in the regression. 
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Table S2: Sensor RMS differences with respect to mouthguard time traces at sensor location 




Linear Acceleration 

Angular Velocity 

Angular Acceleration 

Device 

Parameter 

RMS, g 

NRMS, % 

RMS, rad/s 

NRMS, % 

RMS, rad/s ^ 

NRMS, % 

Skin Patch 

AP Translation/ 
Coronal Rotation 

1.5 (0.3) 

20(5) 

1.2 (0.6) 

78 (29) 

182 (114) 

44 (27) 


LR Translation/ 
Sagittal Rotation 

2.3 (0.6) 

81 (30) 

2.2 (0.6) 

48 (16) 

786 (170) 

120(62) 


IS Translation/ 
Horizontal Rotation 

4.1 (0.8) 

44(7) 

7.5 (2.2) 

700 (380) 

1772 (516) 

860 (500) 


Magnitude 

1.7 (0.4) 

18 (3) 

5.7 (1.9) 

150(76) 

1539 (552) 

290 (230) 

Skull Cap 

AP Translation/ 
Coronal Rotation 

5(2.7) 

47 (40) 

4.8(2.5) 

370 (280) 

1858 (1290) 

460 (340) 


LR Translation/ 
Sagittal Rotation 

5.5 (2.6) 

144 (90) 

4(2.5) 

83 (52) 

1349 (823) 

160(60) 


IS Translation/ 
Horizontal Rotation 

7.3 (2.4) 

79 (20) 

2.6 (1.5) 

330 (450) 

945 (565) 

590(660) 


Magnitude 

6.7 (3.3) 

60 (40) 

4.6 (3.5) 

110 (80) 

2233 (1379) 

500(140) 


Note: average RMS and NRMS errors are reported \A/ith standard deviation in parentheses. 
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Figure S1: Sensor peak measurements projected to head CG. Here we 
show the peak skin patch/skull cap measurements in each individual 
degrees of freedom with respect to the mouthguard reference, in ad¬ 
dition to the magnitude peaks. Note that for individual axes, the sign 
of the peak is taken into account. Peaks in the 1st or 3rd quadrants 
indicate agreement in sign between sensor and reference, while peaks 
in the 2nd or 4th quadrants indicate opposite signs between sensor and 
reference. 
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